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We demonstrate that networks of locally connected processing units with a primitive learning 
capability exhibit behavior that is usually only attributed to quantum systems. We describe 
networks that simulate single-photon beam-splitter and Mach-Zehnder interferometer exper- 
iments on a causal, event-by-event basis and demonstrate that the simulation results are in 
excellent agreement with quantum theory. We also show that this approach can be generalized 
to simulate universal quantum computers. 

KEYWORDS: computer simulation, quantum theory, quantum interference, quantum computation, ma- 
chine learning 
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1. Introduction 

Computer simulation is a powerful methodology to 
model physical phenomena.^ However, some of the most 
fundamental experiments in quantum physics^' have not 
been simulated in the event-by-event manner in which 
the experimental observations are actually recorded.'* In 
experiments the detection of events appears to be ran- 
dom,^' in a sense which, as far as we know, has not 
been studied systematically. Quantum theory gives us a 
recipe to compute the frequency of the observation of 
events but it does not describe individual events, such as 
the arrival of a single electron at a particular position on 
the detection screen.'^' Reconciling the mathematical 
formalism (that does not describe single events) with the 
experimental fact that each observation yields a definite 
outcome is often referred to as the quantum measure- 
ment paradox. This is a central, fundamental problem in 
the foundation of quantum theory.^' Therefore, it is 
not a such a surprise that within the framework of quan- 
tum theory, no algorithm has been found to perform an 
event-based simulation of quantum phenomena. 

From a computational viewpoint, quantum theory pro- 
vides us with a set of rules (algorithms) to compute 
probability distributions.^'^'*" Therefore we may won- 
der what kind of algorithni(s) we need to perform an 
event-based simulation of the experiments^''^ mentioned 
above without using the machinery of quantum the- 
ory. Evidently, the present formulation rules out any 
method based on the solution of the (time-dependent) 
Schrodinger equation and we have to step outside the 
framework that quantum theory provides. 

In this paper we demonstrate that locally-connected 
networks of processing units with a primitive learn- 
ing capability are sufficient to simulate deterministically 
and event-by-event, the single-photon beam splitter and 
Mach-Zehnder interferometer experiments of Grangier et 
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al.^ We also show that this approach can be generalized 
to simulate universal quantum computation by a deter- 
ministic event-by-event process. Thus, the method we 
propose can simulate wave interference phenomena and 
many-body quantum systems using classical, particle-like 
processes only. 

Our results suggest that we may have discovered a 
procedure to simulate quantum phenomena using causal, 
local, deterministic and event-based processes. Our ap- 
proach is not an extension of quantum theory in any 
sense and is not a proposal for another interpretation 
of quantum mechanics. The probability distributions of 
quantum theory are generated by a deterministic, causal 
learning process, and not vice versa.® 

2. Deterministic Learning Machine (DLM)** 

In quantum physics, an event corresponds to the detec- 
tion of a photon, electron, and the like. In our simulation 
approach an event is the arrival of a message at the input 
channel of a processing unit. This processing unit typi- 
cally contains two DLMs (described below). We use the 
diagram of a DLM-based processor that performs the 
event-by-event simulation of single-photon beam split- 
ter, as shown in Fig. 1, to describe the operation of the 
different components of the processor. The applications 
to quantum computations presented later demonstrate 
that the structure of the DLM-based processor is in fact 
generic. 

In Fig. 1, the presence of a message is indicated by an 
arrow on the corresponding line. The first component, 
DLM 1, "learns" about the occurrence of an event on 
one of its two input channels that we label with and 
1. For brevity, we refer to an event on channel (1) as a 
(1) event. The second component transforms the data 
stored in DLM 1 and feeds the result into DLM 2. DLM 
2 "learns" this data. Finally the learning process itself is 
used to determine whether DLM 2 responds to the input 
event by sending out either a or a 1 event. None of 
these components makes use of random numbers, hence 
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Fig. 1. Diagram of the network of two DLMs that performs a deterministic simulation of a single-photon beam splitter (BS) on an 
event-by-event basis. The solid lines represent the input and output channels of the BS. Dashed lines indicate the flow of data within 
the BS. Input channel receives (j/OiJ/l) = (cosi/'o,sini/'o) with probability po- Input channel 1 receives (2/2,2/3) = (cos'!/'l,sini/)i) with 
probability pi = 1 — po- i?(45°) denotes a rotation of a two-dimensional vector by 45° (see Eq. (3)). 



the name deterministic learning machine. 

UsuaUy, DLM-based simulation algorithms contain 
several DLM-based processors that form a network. In 
this paper we only consider networks of processing units 
in which only one message is traveling through the net- 
work at any time. Thus, the network receives an event 
at one of its inputs, processes the event and delivers the 
processed message through one of its output channels. 
After delivering this message the network can accept a 
new input event. 

2. 1 Description of a DLM 

A DLM is a very simple, c;lassical dynamical system 
with a primitive learning capability. This dynamical sys- 
tem consists of a unit vector (such as x = {xo,X\,X2-, xz) 
in DLM 1 of Fig. 1), a rule that specifies how this vector 
changes when an input event is received, and a rule by 
which the DLM determines the type of output event it 
generates as a response to the input event. The initial 
value of the internal vectors is irrelevant. In simulations, 
we usually use random numbers to initialize the internal 
vectors of all the DLMs in the network. 

We now describe the learning process of a DLM in de- 
tail. The basic idea of the learning algorithm is that 
the DLM minimizes the distance between the input vec- 
tor (discussed later) and the internal vector and that this 
minimization is sufficient to construct particle-like pro- 
cesses that mimic quantum phenomena. However, the al- 
gorithm that we describe below cannot be derived from 



the axioms of quantum theory. After many trials and 
failures, we simply discovered that learning algorithms 
of this type can be used to simulate quantum phenom- 
ena. 

First we consider DLM 1 in Fig. 1. The internal state of 
DLM 1 is represented by the vector x = {xq,x\,X2,xz)- 
DLM 1 can accept two different types of input events, but 
only one at a time. Event carries a message represented 
by a two-dimensional unit vector (t/o, Event 1 carries 
a message represented by a two-dimensional unit vector 
(t/2, yz)- Upon receiving an input event, DLM 1 performs 
the following steps: 

• DLM 1 computes eight candidate internal states 

Wi = {+d 1 — + a^x'^,axi,ax2,0iXz), 



W2 = {—Jl — a'^ + a'^XQ,axi,ax2,axs), 



W3 = {axo,+Jl- +a^x\,ax2,axs), 



■W4 = {axo,—Jl — a^-\-a'^x\,ax2,oiXz), 



W5 = {axo,axi,+Jl - + a^xl^axs), 



We = {axo^axi,— J 1 — + 0^X2, oixz), 



W7 = {axQ,axi,Q.X2,+y '\- — oP' + oP'x'^), 
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wg = {axo,axi,ax2,—\J'^ — cP' + cP'x\).{\) 

The parameter < a < 1 controls the leaniiiig pro- 
cess and is discussed in more detail later. The plus 
and minus sign in front of the square roots is in- 
troduced to allow the vector Wj to cover the whole 
eight-dimensional unit sphere. 

• If DLM 1 receives an input event of type 
with message {yoiVi), it constructs a vector x = 
(yO) 2/1) ^2,3:3). If DLM 1 receives an input event of 
type 1 with message {1/2 tUs), it constructs a vector 
X = (xo, xi, 2/2, 2/3)- DLM 1 determines the update 
rule m that minimizes the cost function 

Cj = -wjx, (2) 

that is, Cm < Cj for j = 1, . . . , 8. 

• DLM 1 updates its internal vector by replacing x by 

• DLM 1 generates a new (internal) event by putting 
the values of its internal vector on its four output 
channels. 

• DLM 1 waits for the arrival of the next input event. 

The transformation stage applies an orthogonal trans- 
formation T to X = {xo,Xi,X2, xa). In general, the precise 
form of the transformation T depends on the particular 
function that the processor has to perform. In the ex- 
ample shown in Fig. 1, the orthogonal transformation 
T takes two pairs of elements from x and performs the 
plane rotation 

m = ( '"'t "'"IV (3) 

^ ' \ sm(p cos 9 J ^ ' 

with (p = -k/A. As we show later, this transformation 
implements the single-photon beam splitter. The result 
x' = {x'q,x'i,x'2,x'^) of this transformation is sent to the 
input of DLM 2. Thus, DLM 2 accepts messages in the 
form of a four- dimensional unit vector. DLM 2 updates 
its internal vector z = {zq, Zi, Z2, z^) according to the 
following procedure: 

• DLM 2 performs computes eight candidate internal 
states 

wi = {+Jl - a'^ + a'^z^,azi,az2,az3), 



■W2 = {-Jl - a'^ + a'^Zg,azi,az2,az3), 



W3 = {azo,+Jl - a"^ + a'^zl,az2,az3), 



W4 = {azo,-Jl-a'^ + a'^zf,az2,az3), 



W5 = {azo,azi,+Jl-a'^ +a'^z^,az3), 



We = {azo,azi,-Jl-a'^ + a'^z^,az3), 



W7 = {azo,azi,az2,+Jl - a"^ + a'^z'^), 



Wg = {azo,azi,az2,-Jl - a"^ + a'^z^). (4) 



• DLM 2 determines the update rule m that minimizes 
the cost function 

C, = -wJx', (5) 
that is, C„i < Cj for j = 1, . . . , 8. 

• DLM 2 updates its internal vector by replacing z by 

• DLM 2 generates an output event of type (1) 
if m = 1, . . . , 4 (5, . . . , 8), carrying the message 
iy'o^y'i) = {zo,zi) {{yo,y'i) = {Z2,Z3)). 

• DLM 2 waits for the arrival of the next input event. 

Comparing the algorithms for DLM 1 and DLM 2, we see 

that they are indentical except for part of the second step 
and the fourth step in which the output is generated. 

2. 2 Dynamic behavior of a DLM 

In general, the behavior of a DLM defined by rides 
Ecis. (1) and (2) or Eqs. (4) and (5) is difficult to analyze 
without the use of a computer. However, for a fixed input 
x' = u, it is clear what a DLM will do. It will minimize 
the cost given by Eq. (5) by rotating its internal vector 
z to bring it as close as possible to u. After a number of 
events (depending on the initial value of z, the input u, 
and a), z will be close to u. However, the vector z does 
not converge to a limiting value because the DLM always 
changes its internal vector state by a nonzero amount. 
It is not difficult to see (and supported by simulations, 
results not shown) that once z is close to u, it will keep 
oscillating about u.^^ Below we analyse this behavior in 
more detail using DLM 2 as an example. The dynamics 
of DLM 1 is the same as that of DLM 2. 

Let us denote by no the number of times the DLM 
selects update rule m = 1, 2 (see Eq.(4)). Writing 

= a^zl = {zo + Sf, (6) 

and assuming that <C a < 1, we find that the variable 
Zo changes by an amount 

5^{l-a'){l-zi)/2zo, (7) 

where we have neglected terms of order S^. Similarily, 
if N is the total number of events then N — no is the 
number of times the DLM selects update rules m ^ 1, 2. 
For j ^ 1, 2, Eq.(4) gives 

wlj=a^z^ = {zo + S'f, (8) 

where we have neglected terms of order 6'^. Hence zo 
changes by 

S' ^-{l-a^)zo/2. (9) 

If z oscillates about u then zq also oscillates about 
Up- This implies that the number of times Zo increases 
times the; increment must approximately be equal to the 
number of times Zo decreases times the decrement. In 
other words, we must have noS + {N — no)5' « 0. As 
Zo ~ uo we conclude that no/N « Uq- Applying the same 
reasoning for the cases where the DLM selects update 
rule m = 3, 4 shows that the number of times the DLM 
will apply update rules m = 3, 4 is proportional to Uq + 
ul 
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At this point, there is not yet a relation between the 
dynamics of the DLM and quantum theory. However, 



let us now assume that po = 



- -A 



{pi ^4 + 4) 



is the probability that a quantum system is observed 
to be in the state (1). In quantum theory, we would 
describe this state by a wave function with complex am- 
plitudes zo + izi {z2 + iz:i). Let is now consider a DLM 
that is learning the four values z = {zq, z^, Z2, Zs)- From 
the foregoing discussion, it follows that once the DLM 
has reached the stationary state in which it oscillates 
about z, the rate at which the DLM uses update rules 
m = 1,2, 3,4 (m = 5,6,7,8) corresponds to the proba- 
bility Po (pi) to observe a (1) event in the quantum 
mechanical system. Thus, the DLM generates and 1 
events in a deterministic manner and because it gener- 
ates (1) events if it selected update rule m = 1,2.3,4 
(5,6,7,8), the rate at which these events are generated 
agrees with the corresponding probabilities of quantum 
theory. As the applications presented later demonstrate, 
this correspondence is all that is needed to perform an 
event-by-event simulation of quantum interference and 
many-body quantum phenomena. 

2.3 Stochastic variant 

The sequence of events that is generated by a DLM 
(network) is strictly deterministic but a simple modifi- 
cation turns a DLM into a stochastic learning machine 
(SLM). The term stochastic does not refer to the learn- 
ing process but to the method that is used to select the 
output channel that will carry the outgoing message. 
As explained earlier, in the stationary regime Xq + x\ 
and + x\ (or Zq -\- zf and Z2 + z^) correspond to the 
probabilities of quantum theory. Thus, a comparision of 
for instance, Xq + xf, with a uniform random number 
< r < 1 gives the probability for sending the mes- 
sage over the corresponding output channel. Although 
the learning process of this processor is still determinis- 
tic, in the stationary regime the output events are ran- 
domly distributed over the two possibilities. Of course, 
the frequencies of output events is the same as that of 
the original DLM-network. Replacing DLMs by SLMs in 
a DLM-network changes the order in which messages are 
being processed by the network but leaves the content of 
the messages intact. 

2.4 Generalization 

In the previous discussion, we considered a DLM-based 
processor (see Fig. 1) that accepts two different types of 
events whereby each event carries a message containing 
two real numbers. This is sufficient to simulate quantum 
phenomena such as single-photon interference but if we 
would like to perform event-by-cvent simulations of more 
complicated quantum systems such as quantum comput- 
ers, a generalization is necessary. From the foregoing de- 
scription of the learning rule of a DLM it is obvious how 
this rule may be generalized to handle an arbitrary num- 
ber Ne of different events of messages of arbitrary (but of 
the same) length N„i: Use a vector of NeNm elements to 
represent the internal state of the DLM and, instead of 
eight candidate rules, compare the cost of 2NeNm candi- 
date rules. Clearly, the construction of DLM-based net- 



works is very systematic and straightforward. 

2.5 Summary 

A DLM responds to the input event by choosing from 

all possible alternatives, the internal state that minimizes 
the error between the input and the internal state it- 
self. This deterministic decision process is used to deter- 
mine which type of event will be generated by the DLM. 
The message contains information about the decision the 
DLM took while updating its internal state and, depend- 
ing on the application, also contains other data that the 
DLM can provide. By updating its internal state, the 
DLM "learns" about the input events it receives and by 
generating new events carrying messages, it tells its en- 
vironment about what it has learned. 

3. Single-Photon Beam Splitter 

In quantum theory, the presence of photons in the in- 
put modes or 1 of a beam splitter is represented by the 
complex-valued amplitudes {ao,ai).^'^^'^^^ According to 
quantum theory, the complex- valued amplitudes {bo,bi) 
of the photons in the output channels and 1 of a beam 
splitter arc given by^' 



bo 
bi 



J_ / 1 i 
V2 U 1 



ao 
ai 



(10) 



Writing ao = y^e*'^" and oi = ^1 — poc^'^^ , the prob- 
ability to observe a photon in output channel (1) is 
given by 



|^^|2 ^ 1 + Vpo(l-po)sin(V>o-V>i) ^^^^ 



1 - \/Po(l - Po) sinjipo - ipi) 



(12) 



Here ijjo and ipi represent the phases of the photons. In 
a quantum theoretical description, this phase is propor- 
tional to the length of the optical path that the photons 
have travelled before they enter the beam splitter.^' ^"^^ 

We now show that the DLM-network shown in Fig.l 
behaves as if it is a single-photon beam splitter. This 
network receives events at one of the two input channels. 
There is a one-to-one relation between each input chan- 
nel and the corresponding input mode of the quantum 
mechanical description. Each input event carries infor- 
mation in the form of a two-dimensional unit vector. Ei- 
ther input channel receives {yo,yi) = (cosiAcsin^o) 
or input channel 1 receives (2/2,2/3) = (cosV'ijSin^/ji). In 
terms of the single-photon experiments of Grangier et 
al.,^ an event corresponds to the arrival of a photon at 
channel (1) with phase ipo (V'l) of the beam splitter 
(see Fig. 1). 

The input message is fed into the DLM-network de- 
scribed in Section 2. The purpose of DLM 1 is to trans- 
form the information contained in two-dimensional input 
vectors (of which only one is present for any given input 
event), into a four-dimensional unit vector. The internal 
vector X of DLM 1 learns about the amplitudes (ao, ai): 
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Fig. 2. Simulation results for the beam sjjlittcr sliown in Fig. 1. Markers give the simulation results for the normalized intensity 
No/ {No + Ni) in output channel as a function oi <p = Jpo — Jpi ■ Open circles: po = 1; Bullets: po = 0.5; Open squares: po = 0.25. The 
results of quantum theory |6oP (see Eq. (11)) are represented by lines. 



In the stationary regime we have 

xi w y/posin-tpo, 

X2 « \/l -po COSV-I, 

X3 ~ ^1 -po sinV'i. (13) 

The four-dimensional internal vector of this device is split 
into two groups of two-dimensional vectors {xq,X3) and 
(x2, xi) and each of these two-dimensional vectors is ro- 
tated by 45°. Put differently, the four-dimensional vector 
is rotated once in the (l,4)-plane about 45° and once in 
the (3,2) plane about 45°. The order of the rotations is 
irrelevant. Physically, this transformation corresponds to 
the reflection of the photons by 45° at the beam splitter. 
The resulting four-dimensional vector is then sent to the 
input of DLM 2. The internal vector z of DLM 2 learns 
about the amplitudes {bo,bi): In the stationary regime 
we have 

zo ~ VPocos'^o - - Posinipi, 
zi « y/p^sinipo + -pocosV'i, 
Z2 ~ -PocosV'i - •v/PoCOS'^o, 

Z3 ~ -v/1 -posinV)! + V^cosV'o (14) 

DLM 2 sends {zq, zi) through output channel if it used 
rule m = 1,2,3,4 (see Eq. (5)) to update its internal 
state. Otherwise it sends {z2, Z3) through output channel 
1. 



In Fig. 2 we present results of discrete-event simula- 
tions using the DLM network depicted in Fig. 1. We de- 
note the number of (1) events by Nq (Ni) and the 
total number of events hy N = No + Ni. The correspon- 
dence with the quantum system is clear: the probability 
for a event is given by |6oP « Nq/N, = Re bo/\bo\ 
and y[ = Im 6o/|^o|- The probability for a 1 event is 
|&i|2 « Ni/N, y'^ = Re&i/I&il and y'^ = Im 
Before the simulation starts, the internal vectors of the 
DLMs are given a random value (on the unit sphere). 
Each data point represents 10000 events. All these sim- 
ulations were carried out with a = 0.99. The simulation 
procedure itself consists of four steps: 

(1) Use two uniform random numbers in the range 
[0,360] to generate ?/^o and ^/^i. 

(2) For fixed values of -00 and ■i/'i, generate 10000 in- 
put events. Input channel receives (t/o,yi) = 
(cosV'O)Sin-0o) with probability po • Input channel 1 
receives (2/2,2/3) = (cos V'l, sin V)!) with probability 
Pi = 1 - Po- 

(3) Count the number of output events A^o (-^1) in chan- 
nel (1), see Fig. 1. 

(4) Repeat steps 1 to 3. For each pair {ipo,tpi), store the 
results for Nq (Ni). 

Plotting No/{No + Ni) and |6oP as a function of </) = 
V'o~V'i yields the results shown in Fig. 2. Actually, there 
is no need to use random numbers to generate ipo and 
■01. In Fig. 2, we only used this random process to show 
that the order in which we pick tpo and tpi is irrelevant. 
Random processes enter in the procedure to generate the 




Fig. 3. Diagram of a DLM network that simulates a single-photon Mach-Zchndcr interferometer on an event-by-event basis. -'^-'^ The 
DLM network consists of two BS devices (see Fig. 1) and two passive devices -R{0o) and R{(pi) (see Eq. (.3)) that ijeiform plane 
rotations by ipo and ipi, respectively. The number of events Ni in channel i = 0, ... ,3 is proportional to the probability for finding a 
photon on the corresponding arm of the interferometer. 



input data only. The DLM network processes the events 
sequentially and deterministically. From Fig. 2 it is clear 
that the output of the deterministic DLM-based beam 
splitter reproduces the probability distributions as ob- 
tained from quantum theory. 

4. Mach-Zehnder Interferometer 

In quantum physics, '^'^ single-photon experiments with 
one beam splitter provide direct evidence for the particle- 
like behavior of photons.^' ^ The wave mechanical char- 
acter appears when one performs single-particle inter- 
ference experiments. We now describe a DLM network 
that displays the same interference patterns as those ob- 
served in single-photon Mach-Zchndcr interferometer ex- 
periments.^ The schematic layout of the DLM network 
is shown in Fig. 3. The network described in Section 3 is 
used for the beam splitters. The phase shift is taken care 
of by the devices R{4>o) and R{4>i) (that do not contain 
DLMs) that perform plane rotations by (po and (pi (see 
Eq. (3), respectively. Clearly there is a one-to-one corre- 
spondence between the components of the DLM network 
and the elements of a physical Mach-Zehnder interferom- 



where 



eter 



2, 13 



According to quantum theory, the amplitudes 
(6o, bi) of the photons in the output modes (A^2) and 
1 (N^) of the Mach-Zehnder interferometer are given 

by2.14.15 



(15) 



/ bo ^ 










ax J 



1 











1 i 
i 1 



(16) 



and ao (ai) denotes the amplitude of the photons in 
the input channel (1). Note that in experiments it is 
impossible to simultaneously measure {No/{No + Ni), 
Ni/{No + Ni)) and (N^/iNo + Nx), N^/iNo + Ni)): Pho- 
ton detectors operate by absorbing photons. 

In experiment,^' there are no photons in input chan- 
nel 1 of the first beam splitter, that is ai = 0. From 
Eq. (15), we see that the phase of the wave function de- 
scribing the photons in input channel of the first beam 
splitter is irrelevant. However, as the photons leave the 
first beam splitter the relation between their phases is 
fixed. This relation can be changed through the opti- 
cal path length for reaching the second beam splitter. A 
change of the optical path length in channel (1) results 
in a phase shift by (po {4>i)- If o-i = 0, the output am- 
plitudes of the Mach-Zehnder interferometer are given 

by 



bo 



bi 



ao|ao| 



ao ao 



-l_i(0o-H0i)/2 



sm ■ 



-igi(<^o+0i)/2^og; 



(17) 
\bi? 



from which we see that the probabilities |6oP and 
depend on (f> = (po — (pi only. 

In Fig. 4 we present a representative selection of simu- 
lation results for the Mach-Zehnder interferometer built 
from DLMs. We assume that input channel receives 
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Fig. 1. Simulation results for the DLM-network shown in Fig. 3. Markers give the simulation results for the normalized intensities as 
a function of (j>o. Open squares: No/{No + Ni); Solid squares: N2/{N2 + N3) for 01 = 0; Open circles: N2/(N2 + N3) for 0i = 30°; 
Bullets: N2/{N2+N3) for 0i = 240°; Asterisks: N3/{N2+N3) for 0i = 0; Solid triangles: N3/{N2 + N3) for ^1 = 300°. Lines represent 
the results of quantum theory. 



(2/0,2/1) = (cos -00, sin -00 ) with probability one and that 
input channel 1 receives no events. This corresponds to 
(oo,ai) = (cos'00 + ismipo,0) in Eq. (15). We use uni- 
form random numbers to determine ^q. As in the case 
of the beam splitter, we only use this random process 
to show that the order in which we pick ipo is irrele- 
vant. In all these simulations a = 0.99. The data points 
are the simulation results for the normalized intensity 
Ni/{NQ + N\) for i=0,2,3 as a function oi (j) = (f)Q — (f)\. 
Each data point represents 10000 events (A'o + A^i = 
iVa + iVg = 10000). Initially the rotation angle 0o = 
and after each set of 10000 events, (?io is increased by 
10°. Lines represent the corresponding results of quan- 
tum theory.^" From Fig. 4 it is clear that the determinis- 
tic, event-based DLM network generates events with fre- 
quencies that are in excellent agreement with quantum 
theory. 

5. Universal quantum computation 

It has been shown that an arbitrary unitary operation, 
that is, the time evolution of a quantum system, can be 
written as a sequence of singlc-qubit operations and the 
controlled-NOT (CNOT) operation on two qubits.^''' 
Therefore, in principle, single-qubit operations and the 
CNOT operation arc sufficient to construct a univer- 
sal quantum computer or to simulate any quantum sys- 
tem. In this section, we present results of event-based 
simtilations of single qubit operations and a two-qubit 
quantum circuit containing the CNOT operation to illus- 
trate that DLM-based networks can be used to simulate 
universal quantum computers. 



The state vector of a two-qubit system can be written 

as7.14,17 

1$) = ao|0)i|0)2 + ai|l)i|0)2+a2|0)i|l)2 

+a3|l)i|l)2 

= ao|00) + ai|01) +a2|10) + a3|ll) 

= oo|0)+ai|l)+a2|2)+a3|3), (18) 

where ao , . . . 03 are the amplitudes of the four different 
states and |0)i and |l)j represent the and 1 state of the 
z-th qubit, respectively. For convcnvicncc, in the last line 
of Eq.(18), we represent the basis states of the two-qubit 
system in decimal notation, that is |00) = |0), |01) = |1), 
|10) = |2), and 111) = |3).i7 

5.1 CNOT gate 

By definition, the CNOT gate flips the target qubit if 
the control qubit is in the state |1).^^ If we take qubit 1 
(that is, the least significant bit in the binary notation 
of an integer) as the control qubit, we have 

CNOT|$) = ao|0)i|0)2 + a3|l)i|0)2 + a2|0)i|l)2 

+ai|l)i|l)2 

= ao|00)-Fa3|01)-Fa2|10)-Fai|ll) 

= ao\Q) + as\l)+a2\'l)+ai\2.). (19) 

The schematic diagram of the DLM-network that 
performs the CNOT operation on an event-by-event 
(particle- by-particle) basis is shown in Fig. 5. Concep- 
tually the structure of this network is the same as in 
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Fig. 5. Diagram of a DLM-based processor that simulates a CNOT gate on an event-by-event basis. 




Fig. 6. Quantum circuit representation of two equivalent CNOT operations. The dot and the cross on the line denote the control and 
target qubit, respectively. The square boxes labeled H represent Hadamard gates. 



Table I. Simulation results for the DLM-network shown in Fig. 6, demonstrating that the network reproduces the results of the 
corresponding quantum circuit, that is, a CNOT operation in which qubit 2 is the control qubit and qubit 1 is the target qubit. The 
first half of the events are discarded in the calculation of the frequencies fi for observing an output event of type i = 0, 1, 2, 3. For 200 
events or more, the difference between the event-based simulation results and the corresponding quantum mechanical probabilities is 
less than 1%. 



Processor 


Number of events 


Qubit 1 


Qubit 2 


/o 


/i 


/2 


/3 


Deterministic 


100 








0.98 


0.00 


0.00 


0.02 


Deterministic 


100 


1 





0.20 


0.74 


0.01 


0.04 


Deterministic 


100 





1 


0.16 


0.04 


0.00 


0.80 


Deterministic 


100 


1 


1 


0.16 


0.04 


0.72 


0.08 


Deterministic 


200 








1.00 


0.00 


0.00 


0.00 


Deterministic 


200 


1 





0.00 


1.00 


0.00 


0.00 


Detoriiiiiiistit: 


200 





1 


0.00 


0.00 


0.00 


1.00 


DcLciiiiiiiiblic 


200 


1 


1 


0.01 


0.00 


0.9!) 


0.00 


Stochastic 


2000 








0.965 


0.015 


0.010 


0.010 


Stochastic 


2000 


1 





0.007 


0.970 


0.012 


0.011 


Stochastic 


2000 





1 


0.010 


0.008 


0.016 


0.966 


Stochastic 


2000 


1 


1 


0.005 


0.016 


0.963 


0.016 



the case of the Mach-Zehnder interferometer. As in- 
put to the DLM-network we now have four (0,1,2 or 
3) instead of two different types of events. Each event 
carries a message consisting of two real numbers y = 



{yo,yi), corresponding to the quantum mechanical am- 
plitudes ao, . . . 03. The internal state of each DLM is rep- 
resented by a unit vector of eight real numbers x = 
{xo, ■ ■ ■ ,xr) and there are 16 candidate update rules 
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Table II. Same as in Table I except that the control parameter a = 0.999 instead of a = 0.99 and that ten times as many event were 
generated. The results for the deterininistic simulations are exact within tliree-digit accuracy and have therefore been omitted. 



Processor 


Number of events 


Qubit 1 


Qubit 2 


/() 


/l 


/2 


/3 


Stocliastic 


20000 








0.995 


0.003 


0.001 


0.002 


Stochastic 


20000 


1 





0.002 


0.995 


0.003 


0.001 


Stochastic 


20000 





1 


0.002 


0.001 


0.002 


0.995 


Stochastic 


20000 


1 


1 


0.001 


0.002 


0.997 


0.001 



{{j = 0. ... 7; Sj = ±1}, see Eq.(4)) to choose from. The 
rule that is actually used is determined by minimizing 
the cost function given by Eq.(5). The transformation 

stage is extremely simple: According to Eq.(19), all it 
has to do is swap the two pairs of elements {x2,X3) and 

Instead of presenting results that show that a DLM- 
processor correctly simulates the CNOT operation on an 
event-by-event basis, we consider the more complicated 
network of four Hadamard gates and one CNOT gate 
shown in Fig. 6.^^ Quantum mechanically, this network 
acts as a CNOT gate in which the role of control- and tar- 
get qubit have been interchanged.^^ For the correspond- 
ing DLM-network to work properly it is essential that the 
event-based simulation mimics the quantum interference 
(generated by the Hadamard gates) correctly. 

5.2 Hadamard operation 

The Hadamard operation H is the single-qubit opera- 
tion defined by^^ 



H 



V2 ( 



(20) 



Disregarding phase factors, it performs the same opera- 
tion as a beam splitter. 

The structure of a DLM-processor that performs a gen- 
eral single-qubit operation is identical to the one shown 
in Fig. 1. The only difference is in the transformation 
stage. To implement the Hadamard operation, we use 
the transformation matrix T (see Fig. 1) 



1 

71 



V 



\ 

1 



-1 J 



(21) 



5.3 Simulation results 

In Table I, we present simulation results for the DLM- 
network shown in Fig. 6. Before the first simulation starts 
we use uniform random numbers to initialize the inter- 
nal vectors of the DLMs (ten vectors in total). All these 
simulations were carried out with a = 0.99. From Ta- 
ble I, it is clear that, also for a modest number of events, 
the network reproduces the results of the corresponding 
quantum circuit, that is, a CNOT operation in which 
qubit 2 is the control qubit and qubit 1 is the target 
qubit.i'^ 

As an illustration of the use of SLMs, we replace all 
the DLM 2's by SLMs in the DLM implementation of the 
circuit shown in Fig. 6 and repeat the simulations. Prom 
Tables I and II, we conclude that the randomized version 
generates the correct results but significantly more events 



are needed to achieve similar accuracy as in the fully 
deterministic simulation. 

5.4 Technical note 

All simulations that we presented in this section have 
been performed for a = 0.99. Prom the description of 
the learning process it is clear that a controls the rate 
of learning or, equivalently, the rate at which learned 
information can be forgotten. Furthermore it is evident 
that the difference between a constant input to a DLM 
and the learned value of its internal variable cannot be 
smaller than 1 — a. In other words, a also limits the 
precision with which the internal variable can represent 
a sequence of constant input values. On the other hand, 
the number of events has to balance the rate at which 
the DLM can forget a learned input value. The smaller 
1 — a is, the larger the number of events has to be for 
the DLM to adapt to changes in the input data. 

We use the example of this section to illustrate the ef- 
fect of changing a and the total number of events N. In 
Table II we show the results of repeating the procedure 
used to obtain the data shown in Table I but instc;ad of 
a = 0.99 we used a = 0.999 and adjusted the number of 
events accordingly. As expected, the difference between 
the simulation data and the results of quantum theory 
decreases if 1 — a decreases and the number of events in- 
creases accordingly. Comparing Table I with Table II it 
is clear that the decrease of this difference is roughly pro- 
portional to the inverse of the square root of the number 
of events. 

6. Discussion 

We have proposed a new procedure to construct algo- 
rithms that can be used to simulate quantum processes 
without solving the Schrodinger equation. There is a one- 
to-one correspondence between the components of the 
network and the processing units and the physical parts 
of the experimental setup. Furthermore, only simple ge- 
ometry is used to construct the simulation algorithm. In 
this sense, the simulation approach we propose satisfies 
Einstein's criteria of realism and causality.^ 

An analogy may be helpful to understand the concep- 
tual difference between the conventional description of 
quantum theory and the event-based approach proposed 
in this paper. It is well known that an ensemble of sim- 
ple, symmetric random walks may be approximated by a 
diffusion equation (for vanishing lattice spacing and time 
step). Also here we have two options. If we are interested 
in individual events, we have no other choice than to sim- 
ulate the discrete random walk. However, if we want to 
study the behavior of many random walkers, it is compu- 
tationally much more efficient to solve the corresponding 
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diffusion problem. The latter deseribes the outcome of 
(infinitly) many individual events but does not provide 
information about individual events. The random walk 
is the fundamental mechanism that gives rise to diffusion 
behavior. In this sense, the DLMs described in this paper 
may be regarded as building blocks for a dynamic, deter- 
ministic, local and causal system that generates individ- 
ual events in such a manner that the collective behavior 
of these events is described by quantum theory. 

It may be of interest to compare our approach with 
stochastic wavefunction methods. ^^"^'^ Instead of solv- 
ing the equation of motion of the density matrix, these 
methods solve stochastic differential equations for an en- 
semble of independent realizations of pure states. Typ- 
ically, these methods are used to study open quantum 
systems in which a small number of degrees of freedom 
is coupled to a large reservoir. An attempt to use a vari- 
ant of the stochastic wavefunction method to perform 
an event-by-event study of photon emission was reported 
in Ref.^° In stochastic wavefunction methods, the wave 
function evolves in time according to the time-dependent 
Schrodingcr equation. An uncorrelated random process 
interrupts this evolution to project the wave function 
(that is, make a quantum jump) onto another normal- 
ized state. This evolution of the wave function is similar 
to the change of the internal state of a SLM if we consider 
one isolated event that is processed by the transforma- 
tion stage T and a SLM (see Fig. 1). However, as a SLM 
is capable of learning from previous events, the process 
of generating output events is non-Markovian, this simi- 
larity being very superficial. The fundamental difference 
between the two approaches can also be seen as follows. 
In the stochastic wavefunction method, we can calculate 
the time evolution of each member of the ensemble in 
parallel, at least in principle. In the DLM-approach, this 
is impossible: To exhibit quantum mechanical behavior, 
it is imperative that the DLM-network processes events 
in a sequential manner. In the fully deterministic DLM- 
approac;li (that is, without the randomizing feature of the 
SLM), there is no stochastic process at all. Therefore, 
there also is no relation between the stochastic wave- 
function method and the deterministic, machine-learning 
approach discussed in this paper. 

In conclusion, we have shown that single-particle quan- 
tum interference and quantum computers can be simu- 
lated on an event-by-event basis using local and causal 
processes, without the need of concepts such as wave 
functions or particle-wave duality. 
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